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■ The one-loop corrections to the bottomonium mass spectrum due to the finite 
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, charm mass are evaluated in the framework of the relativistic quark model. The 
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obtained corrections are compared with the results of perturbative QCD. 
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'-^ ' Recently effects of the finite mass of the charm quark in determining the bottomonium 



mass spectrum and the b quark mass were considered within perturbative QCD |2|, ^ 
Here we study these effects in the framework of the relativistic quark model taking into 
account the confining quark interaction. Since the bb bound system is rather far from 
being a Coulombic one (even in the ground state) we expect the size of these effects to 
differ noticeably from that predicted within perturbative QCD, where they were found to 
be substantial 0]. This report may be considered as completing our previous paper 0, 
where the mass spectra of heavy quarkonia were calculated with the account of the one-loop 
radiative corrections. Our model was also successfully applied for describing the heavy-light 
meson mass spectra and radiative and weak decays of heavy mesons 0]. 

In papers ^ a meson is described by the wave function of the bound quark-antiquark 
state, which satisfies the quasipotential equation of the Schrodinger type in the center-of- 
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mass frame: 



^m{p)= J j^,Vip,q-M)^M{q), (1) 



\ 2iJ,ji 2iir) J (27r 

where the relativistic reduced mass is 

M4 - {ml - mlf 

" w ' ^'^> 

and lP'{M) denotes the on-mass-shell relative momentum squared 

Here mi 2 and M are quark masses and a meson mass, respectively. For the bb bound system 
(bottomonium) mi = m2 = m^ and Eqs. (0), (0) take the form 

/i« = — , 62(M) = — -m^. (4) 

The kernel ^(p, q; M) in Eq. (0) is the quasipotential operator of the quark-antiquark 
interaction. It is constructed with the help of the off-mass-shell scattering amplitude, pro- 
jected onto the positive energy states. An important role in this construction is played by 
the Lorentz-structure of the confining quark-antiquark interaction in the meson. In con- 
structing the quasipotential of the quark-antiquark interaction we have assumed that the 
effective interaction is the sum of the usual one-gluon exchange term and the mixture of 
vector and scalar linear confining potentials. The quasipotential is then defined by P] 

V{p, q; M) = Ui{p)u2{-p)V{p, q; M)u,{q)u2i-q) , (5) 

with 

V(p, q; M) = ^asD,Ml^l2 + KL (k)!^^^;^ + Konf (k), (6) 

where as is the QCD coupling constant, D^i, is the gluon propagator in the Coulomb gauge 
and k = p — q; 7^ and u{p) are the Dirac matrices and spinors. The effective long-range 
vector vertex is given by 

r^(k) = 7^ + r = (0,k), (7) 
zm^ 

where k is the Pauli interaction constant characterizing the nonperturbative anomalous chro- 
momagnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 

Konf(r) = (1 - ^)Ar + 5, Konf(0 =eAr, (8) 
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reproducing 



Konf(r)=KonfW+K 



(r) = Ar + B 



(9) 



where e is the mixing coefficient. 

The quasipotential for the heavy quarkonia, including retardation, one-loop radiative 
corrections and expanded in p^/m?, can be found in Ref. All the parameters of our 
model, such as quark masses, parameters of the linear confining potential, mixing coefficient 
e and anomalous chromomagnetic quark moment k, were fixed from the analysis of heavy 
quarkonium spectra |^ and radiative decays [0. The quark masses = 4.88 GeV, rric = 
1.55 GeV, nis = 0.50 GeV, mu,d = 0.33 GeV and the slope of the linear potential A = 0.18 
GeV^ in our model have the fixed values which agree with the usually accepted. The constant 
term B = —0.16 GeV may be adjusted producing the overall level shift. In Ref. p] we have 
considered the expansion of the matrix elements of weak heavy quark currents between 
pseudoscalar and vector meson ground states up to the second order in inverse powers of the 
heavy quark masses. It has been found that the general structure of the leading, ffist, and 
second order l/mq corrections in our relativistic model is in accord with the predictions of 
heavy quark effective theory. The heavy quark symmetry and QCD impose rigid constraints 
on the parameters of the long-range potential in our model. The analysis of the ffist order 
corrections fixes the value of the Pauli interaction constant k = —1 [Q. The same value of k 
was found previously from the fine splitting of heavy quarkonia ^Pj- states IQ . The value of 
the parameter characterizing the mixing of vector and scalar confining potentials, e = —1, 
was found from the analysis of the ^/rriQ corrections |^ and from considering radiative 
decays of heavy quarkonia 0] . 

The perturbative heavy quark potential to two loops including the effects of massive 
loop quarks can be found in Ref. 0. In our model we take into account only one- loop 
contributions since the uncertainty brought by the confining interaction is larger then the 
two-loop contributions. The one-loop correction with the finite c quark mass to the Coulomb 
potential can be included in the form f^: 




r 



ay(r, m) 



TX 
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Ca r31 22 , , ,1 5^ Tp 
wi(r,m,/x) = — |y + y [ln(/xr) +7b]| - -Tp + — 



In — -— + 2Ei[JaQmr) 



,(10) 



where 

7 I rvn / 

-X 



the Euler constant 7^; = 0.5772 and Ca = 3, = 4/3, = 1/2 in QCD. A simple 
representation for the vacuum polarization operator with ao = 5.2 was used 0. Subtracting 
from f i(r, m, /x) its value at m = 0, namely 

31 5 /II 2 \ 
i;i(r, 0, ^^) = —Ca - -T^ + [-^Ca - -Tpj [H^^r) + ^e] , 

we obtain 

2 

At;i(r, m, /i) = Vi{r, m, /i) - t;i(r, 0, /i) = -Tp [ln(v/aj^mr) + 7^; + Ei{^mr)\ . (11) 

Thus the one-loop correction to the static QQ potential in QCD due to the finite c quark 
mass reads as: 

AF(r, m,) = --^i^At;i(r, m„ /^) = [ln(v/^mer) + 7^ + Ei(v/^m,r)] . (12) 

o nr vrr y 

The charm mass correction to the bottomonium mass spectrum is then given by 

AM = {AVmJ , asirub) = 0.22, rUc = 1.55 GeV. (13) 

For averaging we use both the wave functions obtained in calculating the heavy quarkonium 
mass spectra |^ and the Coulomb wave functions. The Coulomb averaging is carried out 
also for as = 0.3 sometimes used in the perturbative QCD description of bottomonium. The 
numerical values of this spin-independent correction for different bb states are presented in 
the Table H 

The Table | shows that for a fixed value of the averagings with and without confining 
potential substantially differ especially for the excited states. For growing n = + L + 1 
the values of (AV^^) slowly decrease in our model whereas for the Coulomb potential they 
fall rapidly. The growth of {SEf,i)^l = (AKnc)coui IB (Table I) is evoked by fast 

increasing values of as{fi)- For these values of as{fi) we also reproduce all the corrections 
(SEf^lY^^^ obtained in |Q. The values of the scale jj were fixed in Ref. from the stability 
conditions containing the corrections 6mb to the b quark mass which are absent in our model. 
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TABLE I: Charm mass corrections to the masses of bottomonium {rir + l}L(n = rir + L + l) states 
(in MeV). 



State 


15(1) 


1P(2) 


25(2) 


1D{3) 


2P(3) 


35(3) 




-12 


-9.3 


-8.7 


-7.6 


-7.5 


-7.2 




-9.5 


-4.2 


-3.8 


-2.3 


-2.2 


-2.1 




-20.7 


-9.7 


-8.8 


-5.5 


-5.2 


-4.9 




-14.3 


-22.1 


-21.9 




-49 


-40.5 




0.277 


0.437 


0.452 




0.733 


0.698 



We estimate the uncertainty in calculating the spin-averaged mass spectrum of the bot- 
tomonium to be few MeV. Some of its sources are the uncertainties in the quark masses 
and confining potential as well as higher order relativistic and radiative corrections. The 
present knowledge of the spin- independent part of the heavy quark potential is rather un- 
certain since it strongly depends on the Lorentz structure of the confining potential. The 
calculated corrections partly may be absorbed in the definition of the constant term in the 
static potential (^). Nevertheless they may become essential in future when the long-range 
QQ interaction will be known better. Then the heavy quark mass should be treated more 
self-consistently not as a fixed phenomenological parameter but as a scale-dependent one 
normalized from the heavy quarkonium mass spectrum. 
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